Search for Lepton Flavor Violation in K + Decays 
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A search for lepton flavor violating decays, 



K+ 



/' 



"e + 7r and tt' 



o 



K 



7r e fj, 



e + /i , was performed using the data collected in E865 at the Brookhaven 



Alternating Gradient Synchrotron. No signal was found in any of the decay modes. At the 90% 
confidence level, the branching ratios are less than 3.0 x 10" 9 , 6.4 x f0" 10 , 5.2 x 10 -10 , 5.0 x 10~ 10 
and 3.4 x 10 -9 respectively. 
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The apparent lepton flavor conservation observed so 
far in particle physics is conveniently accommodated in 
the Standard Model if the neutrino masses are zero. Such 
symmetry can be broken by new physics at a higher en- 
ergy scale, such as Technicolor or Supersymmetry, or by 
neutrinos having Majorana masses. Extensive experi- 
mental efforts have been devoted to searches for lep- 
ton flavor violating kaon decays, K\ — > /i ± e =F Q and 
K + — » 7T + /i + e~ p). In this letter, we report the results 
of a search for — > /i + /i + 7r~ {K^ w ), K + — » e + e + 7r~ 
{K ee „), K+ -> /x+e+TT" (K^) and K+ -» 7r+e+ M - 
(Kire/j.)- Unlike K + — » ir + /.i + e~, which only violates lep- 
ton flavor conservation, these decays also violate genera- 
tion number conservation. In addition, the first three de- 
cays violate total lepton number conservation. and 
K ee7r can proceed by the same mechanism as neutrino- 
less double /3-decays of nuclei if neutrinos have Majorana 
masses. Although the first generation is well explored 
in neutrinoless double /3-decays, provides a unique 

channel to search for effects of Majorana neutrinos in the 
second generation j|. 

The previous searches for K ee7T , Kneit and K„ efJ , were 
performed at CERN 25 years ago @. At the 90% con- 
fidence level (C.L.), the branching ratios were found to 
be Br(K+ -> e+e+Tr" ) < 1 x 1(T 8 , Br{K+ -> 7r+e+/i~ 
) < 7 x 10~ 9 and Br(K+ -> n+e+ir- ) < 7 x 1CT 9 . In 
a reanalysis of data of a 1968 bubble chamber experi- 
the best limit on K ua7: was determined to be 



ment 

Br{K~+ -> fi+n+n- ) < 1.5 x 10~ 4 at the 90% C.L. |§. 

Experiment E865 at the Brookhaven Alternating Gra- 
dient Synchrotron was primarily designed to search for 
K + — > n + n + e~ [||. Because of its excellent capability 
in kinematic reconstruction and particle identification of 



K + decays to three charged particles, it has been ex- 
ploited to study other decays such as K + — ► 7r + e + e~ 
§, K+ -> e+ue+e- and K+ -> fi+ve+er. In 1997, 
two special data sets were collected to study K + — > 
ir + /j+ fj,~ (K^y) and K+ K+ir-e+v (A e4 ). Over 400 
K^nn events Q and 400,000 K e ± events were observed. 
We use the former to search for A" MM7r , and the latter for 
A ee7r , A^-e^ and A^ e7r . 
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The detector (Fig. and its performance has been 
described in other publications ^|J^-|^]. The apparatus 
resided in an unseparated 6 GeV beam directly down- 
stream of a 5m long evacuated decay volume. The 
charged particles from K + decays were first separated 
in charge by a dipole magnet, then momentum ana- 
lyzed in a spectrometer system consisting of proportional 
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chambers(Pl-P4) and another dipole magnet. Particle 
identification was achieved by two sets of Cerenkov coun- 
ters (CI) upstream and (C2) downstream of the spec- 
trometer magnet, and a shashlyk-style electromagnetic 
calorimeter downstream of the spectrometer system, fol- 
lowed by a muon range stack consisting of steel plates 
interleaved with proportional tubes. For this study, the 
Cerenkov counters were filled with methane gas for high 
e ± identification efficiency. 

The trigger hodoscopes were located directly down- 
stream of the first proportional chamber PI (D-hod), 
upstream of the calorimeter (A-hod), and in the mid- 
dle (B-hod) and at the end of the muon stack (C-hod). 
The first level trigger was constructed by requiring two 
charged particles on the right, and one charged particle 
on the left in the A and D hodoscopes and the corre- 
sponding calorimeter modules. In the next trigger level, 
particle identification information was applied. 

The trigger designed for K e 4 accepted events with e + 
but not accompanied by an e~ . Cerenkov light signals 
were required on the right side of both CI and C2, and 
both Cerenkov counters on the left were required to have 
signal below one photo-electron, to suppress events with 
an e~ from the 7r° — > e + e _ 7 decay (Dalitz). 

The trigger designed for decay required one 

muon on each side of the detector. Each muon, for trigger 
purposes, was identified as a spatially correlated coinci- 
dence between the B and C hodoscope hits. 

In the off-line reconstruction, events are required to 
have three charged tracks from a common decay vertex 
in the decay volume, a reconstructed kaon momentum 
consistent with the beam phase space distribution, and 
a timing spread between the tracks consistent with the 
resolution, typically about 0.5ns. Similar to the analy- 
sis of the K vfl ^ events JjJ , a joint likelihood function is 
constructed based on the vertex quality, the kaon mo- 
mentum vector, and the track % 2 . This is used to select 
events with high kinematic quality. 

For events, muons are required to have momenta 

greater than 1.3 GeV/c, go through the muon stack and 
have corresponding hits in B-hod and C-hod. There 
should be sufficient muon chamber hits associated with 
the track, and energy deposition in the shower calorime- 
ter should be consistent with minimum ionizing particles. 
The trigger requirement that there be one muon on the 
left and one on the right is not efficient for this decay 
because positively charged particles tend to populate the 
right side of the detector. The majority of the 
events which would be accepted by the trigger would 
have two fi +, s on the right side of the spectrometer sys- 
tem, and one of the /i + 's crossing to the left in the muon 
system downstream of the calorimeter (see Fig. |l|). In a 
smaller fraction of events, one of the ^ + 's stays on the 
left side throughout the detector, and the other and 
the pion on the right. Monte Carlo simulation shows that 
the trigger acceptance for is a factor of 2.7 smaller 



than that for K^^. 

The background for comes from K + — > 

7r + 7r + 7r~ (K T ), with both 7r + 's misidentified as /i + 's. Al- 
though most of these background events have the recon- 
structed /j^jltt mass much lower than Mk because of the 
mass difference between muon and pion, events with pion 
decays in the spectrometer magnet can result in fifxn 
mass in the signal region. Because those events tend to 
have worse kinematic characteristics, a tight cut on the 
joint likelihood helps to reduce background. 
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FIG. 2. The watt invariant mass distribution for 
candidates. The points with error bars are data, the solid 
line is a fit to an empirical function jjj, and the dashed line 
is a fit that includes a signal at the 90% C.L. upper limit. 

Figure ^ shows the reconstructed fifnr mass distribu- 
tion after a cut on the joint likelihood function. The 
background within the signal region (0.4875 GeV < 
M^^-n < 0.5025 GeV) is estimated by fitting the spec- 
trum with an empirical function used in the K^pp analy- 
sis with the signal region excluded. There are 5 events 
in the signal region where 5.3 background events are ex- 
pected. Using the frequentist approach ||, the upper 
limit on the number of signal events is 4.8 at the 90% 
C.L.. Normalizing to K T , we obtain an upper limit on 
the branching ratio: 

Br{K+ -> h+h+tt-) < 3.0 x 10" 9 (90%C.L.). (1) 

For K + — > e+Tr^/iT events, an e + is required on the 
right side with Cerenkov light associated with the track 
in both CI and C2, and an E/p ratio of at least 0.8. The 
charged pion is required to have no significant signals 
in the Cerenkov counters associated with the track, and 
calorimeter responses consistent with minimum ionizing 
particles or hadronic showers. The /z _ (/i + ) is required to 
be on the left(right), to reach the B-hod, and to have a 
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range in the muon stack consistent with its momentum. 
The minimum momentum for the muon is 0.75 GeV/c. 

The main sources of background for K + — > e + 7r ± /i T 
decays are K e 4, when one of the charged pions is misiden- 
tified as muon, and K T , when one 7r + is mistaken for a 
muon and the other n + misidcntificd as e + . The prob- 
ability of misidentifying a 7r + as fi + is 5% due to pion 
decays and punchthrough. The probability of misiden- 
tifying a 7r + as an e + is l.OxlO -4 . This happens when 
the pion deposits most of its energy in the calorimeter, 
and at the same time there are photoelectrons associ- 
ated with the track, either originating from scintillation 
or random activity. Since the threshold of the Cerenkov 
counters is 3.5 GeV for muons, the high energy muons 
in the beam halo can produce Cerenkov light. To reduce 
this misidentification probability, events with additional 
tracks on the right side, either electrons or high energy 
muons, are rejected from this sample. 

In Fig. ||, data are compared to the Monte Carlo simu- 
lation of the background events from K e 4 and K Tl before 
a tight cut on the joint likelihood function is imposed. As 
can be seen, these two decay modes successfully account 
for the observed background. 
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FIG. 4. Scatter plot of M-khh and joint likelihood func- 
tion for K^e^ candidate events. The box indicates the signal 
region. 
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FIG. 5. Similar to Fig. g but for K„ e 



FIG. 3. 7re/x invariant mass distribution for K ne ^ candi- 
dates before a tight cut on the joint likelihood function. The 
points are data, the dashed histogram is the K e 4 Monte Carlo 
simulation, the dotted histogram is the K T Monte Carlo sim- 
ulation, and the solid histogram is the sum. 

Because of the undetected neutrino in if e 4, the K e i 
background is greatly reduced by requiring the candi- 
dates to have a reconstructed kaon momentum vector 
within the beam phase space. Due to the large difference 
in rest masses, the K T background has a lower recon- 
structed invariant mass. 



Another potential background is K + — > 7r + 7r°, fol- 
lowed by 7T° Dalitz decay. Since the misidentification 
probability of e~ as pT is negligible, it does not con- 
tribute to K^eiL background. For events, the re- 
duction of the remaining Dalitz background to a negligi- 
ble level is achieved by requiring M ee > 50MeV, where 
M ee is the invariant mass of e + and tt~ with the 7r~ mass 
assigned to be the electron mass. 

Figures [| and || are the scatter plots of the invariant 
mass of the reconstructed candidate events vs. the joint 
likelihood function. The boxes indicate the signal region, 
which covers ±3cr in mass, and 80% acceptance in joint 
likelihood function. No signal events are observed. 
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The search for -R" ee7r applies particle identification con- 
ditions as described above. The background comes from 
K e 4 where the tt + is misidentified as an e + , and from 
K T where both 7r + 's are misidentified as e + 's. Because 
of the more significant mass difference between 7r's and 
e's these backgrounds are far away from the K ee7T signal 
region. Figure ^| shows the scatter plot of M eep vs. the 
joint likelihood function. Again, there are no events in 
the signal box. The background events in this plot are 
correctly accounted for by K T and K e4 . 
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Br{K+ -» e+e+Tr") < 6.4 x 10~ 10 
Br(K+ -> 7: + e + fi-) < 5.2 x 10~ 10 
Br(K+ -> /i+e+Tr") < 5.0 x 1CT 10 
B^tt -> e + /z~) < 3.4 x 1(T 9 
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AccKei^AcCs) is the acceptance to K e4 (signal) decay, 
and N s = 2.44, Br Kei = 3.91 x 1CT 5 , N Ke4 = 378,000. 
For 7r efl , the branching ratio of 0.21 is taken into 
account. 

The limits on K^ eyL , K^ elT and K ee7r represent an im- 
provement of more than a factor of 10 over the previ- 
ous searches Q. The upper limit on Br(7r° — > e + /i 
and our result of Br(7r° — > /i + e~) < 3.8 x 10~ 10 
can be compared to the previous best limit of [Br(TT° 
H + e~) + Br{7T° -> ^-e+)\ < 1.72 x 10~ 8 @. The new 
upper limit on K^^Q^q. |l|) is a factor of 50000 better 
than the previous experimental bound. The implications 
of this result are discussed in jll| . 
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FIG. 6. Similar to Fig. W, but for K e 



The observation of no signal event for K + — » 7r + e + /i~ 
implies also a null result for the search of tt° — » 
e + /i~(7r e „), through the decay K + — > 7r + 7r°(/^ 2 )- 

Table | lists the acceptances for the decays of inter- 
est. Normalized to K e4 decay, the null results of these 
searches are expressed in term of the 90% C.L. upper 
limit of the branching ratios, 



Decay 


Acceptance (%) 


7r + 7r~ e + v 


3.93 


e + e + 7r~ 


1.54 


7r + e + fi~ 


1.90 


H + e + -K~ 


1.97 


7T + 7r°, 7T° — > e + fl~ 


1.38 


7r + 7r + 7r - 


6.25 




0.71 



TABLE I. The acceptances for K e4 , K^^, K ee7T , K^en, 
Kt, Kpnn and K 7T 2+'^e t i decays, using Monte Carlo simu- 
lation. Detector efficiencies and effects of all kinematic and 
particle identification cuts are included. 
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